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5.2 Kelvin- Helmholz Instability for continuous shear 
and stratification 

5.2.1 Heuristic reasoning 

Due to viscosity, shear flow exists along the boundary of a jet, a wake or a plume . On the 
interface of salt and fresh water, density stratification further comes into play. When will 
dynamic instability occur? 



Refering to Figure 5.2.1, Consider two fluid parcels, each of unit volume, at levels z and 
z + dz. Let their positions be interchanged. To overcome gravity, the force needed to lift the 
heavier fluid parcel by r\ is 




Figure 5.2.1: Exchanging fluid parcels in a stratified shear flow 
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where Boussinesq approximation is used. After the exchange, the parcels mix with the 
surrounding fluid and attain the average velocity 

(U + U + dU)/2 = U + dU/2 

Therefore the total kinetic energy is 

p{U + dU/2) 2 

The available kinetic energy is the difference between the kinetic energies before and after 
the exchange. 

| {u 2 + (U + dUf - 2(U + dU/2) 2 } = ^dU 2 . 

If the net available kinetic energy exceeds the work needed for the exchange, the disturbance 
will grow and the flow will become unstable, i.e., 

P dU2 ^ A-A 

— - — > —gdpdz 
Let the Richardson number be defined by 



_gdp 

Ri = — ^ (5.2.1) 



dU 

dz 



7 > Ri = V^2 (5-2.2) 



Instabilty occurs if 

4 ""• (f ) 

(Chandrasekar, 1961 ). 

Remark : A slightly more accurate estimate can be made without Boussinesq approxima- 
tion. Before the exchange, the total kinetic energy is 



l -{pU 2 + {p + dp){U + dU) 2 }. 
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After the exchange, the parcels mix with the surrounding fluid and attain the average velocity 

(U + U + dU)/2 = U + dU/2 
but their densities are preserved. Therefore the total kinetic energy is 

^p + p + dp)(U + dU/2) 2 

The available kinetic energy is the difference between the kinetic energies before and after 
the exchange. 

~^dU 2 - UdUdp + ^dpdU 2 
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Ignoring the last term, the necessary condition for instability is 

~^dU 2 - UdUdp + ^dpdU 2 > -gdpdz 

or 

1 -I 7 

1 p dz _|_ J- ^ p dz 



4 4 p r^ 2 



dz / 

On the left-hand side, the third term is negligible compared to the first. The ratioi of the 
second term on the left to the term on the right is 

UdU U*_ 
g dz gL 

where L is the length scale of stratification. As long as the last ratio is very small, the 
criterion Ri < 1/4 still holds. 

Let us confirm the heuristic result but the linearize theory. 



5.2.2 Linearized instability theory for continuous shear and strat- 
ification. 

Let the total flow field be {U + u, w, P+p, p + p) where U, P, p represent the backgraound flow 
(u, w,p, p) the dynamical perturbations of infinitesimal magnitude. The linearized governing 
equations are: continuity: 

u x + Wz = Q (5.2.3) 



incompressiblity: 
where 

and momentum conservation: 



Pt + Up x + wjf = (5.2.4) 



dz 



p (u t + Uu x + wU z ) = -p x (5.2.5) 

p(w t + Uw x ) = -p z - pg. (5.2.6) 

where p denotes the perturbation of density from p. 

Let us follow Miles and introduce a new unknown r\ by enoting p = —p'r), then Eqn. 
(5.2.4) gives 

Vt + Ur] x = w (5.2.7) 

Consider 

r] = F(z)e lk{x - Ct \ (5.2.8) 
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where 

c = uj/k = c r + ici. 
For fixed k the flow is unstable if q > 0, since 

Let 

= {u(z),w{z),P(z),-p'F(z)}e^-^ (5.2.9) 

We get from Eqn. (5.2.7) 

w = ik(U-c)F,, (5.2.10) 

from Eqn. (5.2.3) 

u= -[(17- c)F]', (5.2.11) 

and from Eqn. (5.2.5) 

p (iJfe(l/ - c)u + C/'[ifc(C/ - c)F}) = pik 

or 

P[(U - c)(-)[(U - c)F\ + U\U - c)F] = p, 

hence 

p = p{U-c) 2 F'. (5.2.12) 
Substituting Eqns. (5.2.9), (5.2.10), (5.2.11) and (5.2.12) into Eqn. (5.2.6), we get 

[p(U - c) 2 F']' + p[N 2 - k 2 (U - cf] F = 0, (5.2.13) 
where N is the Brunt-Vaisala frequency defined by 

N 2 = - 9 -^. (5.2.14) 

p dz 

Let the top and bottom be rigid walls, then w — 0. Hence, 

^ = i.e., F = 0, z = 0,d. (5.2.15) 

The argument is unchanged if the top and bottom are at z = oo and z = —oo. Equations 
(5.2.13) and (5.2.15) consititute an eigenvalue problem where c = c r + ici is the eigenvalue. 
If q > 0, instability occurs. 



5.2.3 A necessary condition for instability (J.W. Miles, L. N. 
Howard). 



For brevity we set W = U — c. Miles further introduce G = yWF, so that Eqn. (5.2.13) 
becomes 

1 G = 0. (5.2.16) 



(pWG')' - 
The boundary conditions are 

G(0) = G(d) = 0. (5.2.17) 
Multiplying Eqn. (5.2.16) by G* and integrating by parts 

f* [pW (| G[ | 2 +k 2 | G 1 | 2 ) + \ (pU')' \G\ 2 + p (if/' 2 - iV 2 ) W | 2. |2} <fe = o. 

(5.2.18) 

We now seek the necessary condition for instability, i.e., q 7^ 0. Writing 

W = (U - Cr) - ici W* = (U- c r ) + ia 
and substituting these in (5.2.18), we get 

j*{p(U-c r -ic i )(\G'\ 2 +k 2 \G\ 2 ) 

+\ (PU')' \G\ 2 + P (\u' 2 - N 2 ) (U-c r + i Ci ) I ^ | 2 } dz = 0. 
Separating the imaginary part, we get, if q 7^ 0, 

jf P ((I G' | 2 +A; 2 I G | 2 ) <fe + j% (g(5 - \{U'f) \ ^ | 2 dz = 0. 
For this to be true it is necessary that N 2 < \{U') 2 or 

N 2 1 

This confirms the heurisic result as the necessary (but not sufficient) condition for instability 
(J.W. Miles, L. N. Howard). 



